A069570 


r\V~) — ■ 

Tahte(n}|r78iii 


ASSESSMENT  OF  A JtBTHOD  FOR  ^PREDICTING  THE  DRAG  AND  BOUNDARY-^ YER 
DEVELOPMENT  ON  A TODY  OF  REVOLUTION  AT  ^RO  INCIDENCE 


A «i /atkins 


IfTdRuH 


Summary 

A comparison  between  experimental  data  and  a boundary-layer  and  drag 
prediction  method  for  bodies  of  revolution  at  zero  incidence  and  high  Reynolds 
number  is  discussed.  The  only  physical  input  to  the  prediction  method  is  that 
of  transition  position.  In  general  drag  predictions  agree  well  with  measure- 
ments both  in  trend  and  in  absolute  magnitude  and  predicted  velocity  profiles 
on  afterbodies  with  an  included  semi-angle  of  less  than  about  15°  agree  well  with 
measurements.  A modification  to  the  method  has  been  developed  which  gives  a 
better  agreement  between  predicted  and  measured  velocity  profiles  on  fuller 
afterbodies  provided  that  no  flow  separation  takes  place..  (U) 
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LIST  OF  SYMBOLS 


Profile  drag  coefficient  baaed  on  wetted  area. 

Profile  drag  coefficient  baaed  on  maximum  croas-aectional  area. 
Skin-friction  drag  coefficient  baaed  on  wetted  area. 

Local  akin-friction  coefficient. 

pyx 

Profile  drag  coefficient  baaed  on  (volume V*  * , 

Body  diameter. 

Shape  factor. 

Body  length. 

Length  of  noae  section. 

Length  of  parallel  centrebody. 

Length  of  tail  aection. 

Power-law  index. 

Wake  defect  parameter. 

Body  radiua. 

Reynolds  number  baaed  on  body  length  and  free-atream  velocity. 
Local  velocity. 

Total  velocity. 

Free-atream  velocity. 

Velocity  at  edge  of  boundary  layer. 

Hypothetical  velocity  uaed  in  aection  2.4. 

Distance  from  noae  along  the  axia. 

Distance  from  hull  along  a normal  to  the  hull. 

Angle  of  tangent  to  body  surface  with  respect  to  the  axia. 
Displacement  area. 

Boundary  layer  thickness. 

Displacement  thickness. 

Momentum  area. 


INTRODUCTION 


In  the  development  end  feasibility  studies  of  underwater  vehicles,  methods 
for  predicting  the  drag  and  boundary  layer  development  of  a body  of  revolution 
are  required.  They  can  be  used  to  provide  initial  estimates  of  drag,  possibly 
reducing  the  requirement  for  costly  wind-tunnel  or  towing- tank  tests. 

A widely  used  boundary  layer  and  drag  prediction  method  in  the  UK  has  been 
that  of  Young  £lj  which  was  used  to  produce  the  Engineering  Sciences  Data  Unit 
(ESDU)  data  items  f 2 J . The  accuracy  of  this  method  is  limited  inasmuch  as  it 
only  takes  into  account  Reynolds  number,  transition  position  and  body  fineness 
ratio.  It  does  not  distinguish  between  bodies  with  the  same  fineness  ratio  but 
different  body  contours. 

There  are  now  available  a number  of  methods  for  predicting  the  drag  and 
boundary  layer  development  of  axisyametric  bodies.  The  method  used  in  this 
Report  was  obtained  from  Dr  Myring  of  Salford  University  and  details  of  how  to 
use  it  appear  in  an  earlier  Note  . In  Reference  3,  predicted  boundary  layer 
integral  parameters  on  the  afterbody  are  compared  with  estimates  obtained  from 
measurements.  In  general,  there  is  agreement  between  theory  and  experiment  in 
trend,  if  not  always  in  absolute  magnitude.  In  particular,  on  blunt  afterbodies, 
the  shape  parameter  is  under-predicted  near  the  tip  of  the  body. 

One  of  the  applications  of  the  prediction  method  is  to  determine  scaling 
effects  on  boundary  layer  velocity  profiles.  In  this  case  values  of  the  velocity 
components  are  required  rather  than  integrals  of  them.  In  the  present  Report, 
predicted  and  measured  velocity  profiles  on  slender  afterbodies  agree  well  but 
the  agreement  is  less  favourable  on  afterbodies  with  an  included  semi-angle 
greater  than  about  15° • This  disparity  is  discussed  and  a method  of  overcoming 
it  is  examined  by  comparing  predicted  velocity  profiles  and  boundary  layer 
integral  parameters  with  experiment. 


2.  THEORY  OF  THE  PREDICTION  METHOD 


2.1  Previous  methods. 


All  theoretical  methods  of  drag  prediction  involve  calculation  of  the 
boundary  layer  development  over  the  body  surface.  Young  fl]  showed 
theoretically  that  the  profile  drag  is  proportional  to  the  momentum  area 
of  the  wake  far  downstream  of  the  tail  and  he  related  the  values  at  infinity 
to  those  at  the  tip  of  the  body.  In  the  turbulent  boundary  layer  calculation 
an  experimental  pressure  distribution  was  used  as  input  and  the  shape 
parameter  was  assumed  to  be  constant.  Granville  fV]  used  a one-parameter 
integral  method  which  allowed  the  shape  parameter  to  vary  and  incorporated  a 
criterion  fcr  predicting  transition.  More  recently  Cebici  et  al  [5,  6J 
used  a finite  difference  method  with  a simple  eddy  viscosity  approach  for 
modelling,  the  • Reynolds  stresses.  Hess  & James  [7]  using  this  method 
demonstrated  that  great  care  was  needed  in  extrapolating  from  an  axial 
station  at  or  near  the  tip  to  infinity  in  order  to  obtain  the  correct  value 
of  drag. 

Previous  calculation  methods  suffer  from  the  disadvantage  that  the 
pressure  distribution  over  the  body  surface  needs  to  be  calculated  or 
measured  separately.  There  are  accurate  potential  flow  calculation  methods 


for  doing  this  even  on  sharply  curved  or  pointed  bodies  f&J,  but  these 
do  not  take  the  wake  into  account.  The  usual  technique  in  this 
circumstance  is  to  extrapolate  the  potential  flow  pressure  distribution 
into  the  wake.  A further  drawback  at  the  tail  end  of  the  body  is  that 
the  boundary  layer  thickness  becomes  large  compared  with  the  body  radius 
and  the  usual  assumptions  associated  with  boundary  layer  theory  cease  to 
apply.  Also  on  full  afterbodies,  the  pressure  distribution  predicted 
using  potential  flow  theory  does  not  correspond  to  that  which  the 
boundary  layer  sees.  Nakayama  & Patel  T9 7 overcame  this  by  using  an 
integral  method  adapted  to  calculate  thick  axiisymmeatrcuE  turbulent 
boundary  layers.  After  each  step  in  the  boundary  layer  calculation  over 
the  last  ten  per  cent  of  the  body  length,  a value  of  drag  was  computed 
by  extrapolation  to  infinity  and  the  final  value  chosen  was  the  maximum 
of  those  obtained.  Kerney  & White  fioj  used  a two-parameter  integral 
method  of  Granville  fll J applicable  to  thick  or  thin  boundary  layers. 

2.2  The  present  method. 

A method  which  attempts  to  overcome  most  of  the  previously  mentioned 
disadvantages  has  been  developed  by  Myring  £127*  The  method  couples  the 
boundary  layer  and  potential  flow  calculations  in  an  iterative  cycle  and 
the  potential  flow  calculation  is  carried  out  over  the  body  surface  plus! 
displacement  thickness.  The  entire  calculation  is  extended  into  the  wake 
where  at  some  prescribed  location  the  displacement  surface  is  truncated 
and  suitable  end  conditions  are  imposed.  Apart  from  the  details  of  body 
shape  the  only  inputs  required  are  those  of  transition  position  and 
free-stream  conditions. 

The  potential  flow  calculations  uses  the  well  established  surface  source 
distribution  method  of  Hess  & Smith  £ 8],  The  laminar  part  of  the  boundary 
layer  is  calculated  using  the  method  of  Luxton  & Young  £137  and  the 
turbulent  part  using  Head's  entrainment  method  £l4].  The  far  wake  is 
treated  by  using  an  alternative  entrainment  hypothesis  and  the  near  wake 
by  interpolation  between  the  attached  flow  and  the  far  wake.  The  profile 
drag  is  computed  from  the  value  of  momentum  area  in  the  far  wake  which  is 
extrapolated  to  infinity. 

The  method  also  calculates  the  skin-friction  drag  which  is  obtained 
by  integrating  the  local  skin-friction  over  the  body  surface.  Two  skin- 
friction  laws  are  currently  programmed  into  the  prediction  method.  Myring 
used  the  empirical  relationship  of  Ludwieg  & Tillman  £l5j.  Chappell  [1$1 
pointed  out  that  for  Reynolds  numbers  R,  (based  on  body  length)  above  10' 
the  skin-friction  drag  of  a flat  plate  predicted  using  this  relationship 
begins  to  diverge  from  experimental  data,  and  that  for  slender  bodies  of 
revolution  at  high  Reynolds  number  (R.>108)  the  predicted  total  body  drag 
is  less  than  that  of  a flat  plate  witn  the  same  surface  area  and  transition 
position.  Chappell  £1 7]  modified  the  Ludwieg-Tillman  expression  to  provide 
a better  fit  to  experimental  flat  plate  data  over  a wider  range  of  Reynolds 
number.  The  modified  law  agrees  well  with  flat  jjliteidataeiqcfc^fi  range 

1o5<rl<io9. 


- 8 - 


2.3  Determination  of  velocity  profiles. 

The  method  predicts  directly,  the  momentum  area  0 defined  by 

«-r  (U/U^  (1  - U/U^  (r  + y cos  )dy  (l) 

and  the  shape  parameter  H defined  by 


H = A*/0 

where  a *-r  (l  - U/U^ ) (r  ♦ y coseOdy 


(2) 

(3) 


In  (l)  and  (3)  r denotes  the  body  radius  and  o(  the  angle  of  the  body 
surface  with  respect  to  the  axis.  The  integrations  are  performed  along 
a line  normal  to  the  hull.  In  order  to  obtain  values  of  the  velocity 
component  parallel  to  the  hull,  it  is  assumed  that  the  velocity 
distribution  across  the  attached  boundary  layer  can  be  expressed  in  a 
power-law  form 

o/u,  = (y/£)1/n  <*0 


while  in  the  wake 

U/U  = 1-P 

1 


+ cos 


(5) 


where  (1-2P)  is  the  value  of  U/U  on  the  wake  centre  line.  In  the 
calculation  procedure,  the  parameters  $ (boundary  layer  thickness),  n 
(power-law  index)  and  P (wake  defect  parameter)  are  estimated  from  the 
external  conditions  and  the  boundary-layer  integral  parameters. 


2 . 4 Full  after-body  modification..  . 

As  a result  of  a preliminary  assessment  of  the  method  £*3? , a 
modification  has  been  developed  for  full  after-bodies  following  the  work 
of  Patel  Z l8jl»  Patel  used  a method  requiring  knowledge  of  the  longitudinal 
pressure  gradient  both  at  the  edge  of  the  boundary  layer  and  at  the  wall, 
and  showed  using  hot  wire  measurements  that  they  were  not  the  same.  The 
longitudinal  pressure  gradient  at  the  edge  of  the  boundary  layer  is  related 
to  the  corresponding  velocity.  Let  U denote  the  measured  velocity  at  the 
edge  of  boundary  layer  and  TJ  the  corresponding  hypothetical  velocity 
deduced  from  the  pressure  gradient  at  the  wall.  The  prediction  method  was 
ipodified  by  assuming  that  U^  and  are  related  as  follows: 

til  = U1  £l  - (15  - I*  - l|  )/225|  (0.85<Xi(L<i.i5) 

= U^  elsewhere  (6) 

This  empirical  law  is  based  on  the  measurements  of  Patel  et  al  £l£7* 

Suppose  U.  is  the  external  velocity  predicted  using  the  potential  flow 

calculation.  The  TT  was  used  in  the  boundary  layer  computation  rather 

than  U„ . ' 

1 
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The  location  of  the  natural  transition  position  is  difficult  to 
predict.  Several  different  criteria  have  been  proposed  for  bodies  of 
revolution  and  a review  of  these  was  carried  out  by  Granville  £203  and 
by  Kaups  £2lJ.  According  to  both  authors  none  of  the  available  prediction 
methods  agrees  well  with  experiment  over  a wide  range  of  conditions  and 
body  shapes.  Granville  however  deduced  an  empirical  correlation  of 
transition  measurements  on  streamline  bodies.  Nakayama  and  Patel  [9j 
incorporated  four  different  criteria  in  a boundary  layer  prediction 
method  and  these  gave  results  which  often  agreed  among  themselves,  but 
differed  widely  from  the  observed  values.  The  available  experimental 
evidence  (surveyed  by  McCarthy  et  al  £"223)  suggests  that  the  position  of 
natural  transition  varies  with  nose  shape,  Reynolds  number,  the  level  of 
turbulence  in  the  free-stream  and  the  roughness  of  the  body  surface. 
However,  in  practice,  at  high  Reynolds  number  and  high  prismatic 
coefficient  (4V/7TD2l),  transition  is  likely  to  occur  well  forward  and 
will  normally  be  controlled  by  a tripping  mechanism  in  model  tests. 

For  these  reasons  transition  position  is  prescribed  explicitly  rather 
than  predicted  in  the  present  method. 

2.6  Limitations  of  the  method. 

There  are  situations  where  the  prediction  method  will  not  work 
owing  to  numerical  instabilities  arising  though  the  viscous -•inviscid 
coupling.  The  instabilities  are  caused  by  the  laminar  boundary  layer 
method  becoming  unstable  in  large  favourable  or  adverse  pressure 
gradients.  The  most  significant  practical  example  of  this  is  a body 
with  a flat  nose,  often  used  for  torpedoes.  One  way  to  obtain 
predictions  for  such  a body  is  to  fit  a hemispherical  or  elliptic  nose, 
which  gives  rise  to  added  uncertainties  in  the  accuracy  of  the  results, 
which  are  not  large  if  the  flat  nose  shape  does  not  give  rise  to  laminar 
separation  on  the  nose. 

An  assumption  inherent  in  the  method  is  that  the  flow  does  not 
separate  before  the  tip  of  the  body.  This  means  that  very  blunt  or 
rounded  tails  are  excluded.  Experimental  evidence  mentioned  below 
suggests  that  laminar  separation  from  the  nose  will  only  occur  at  low 
Reynolds  numbers  (typical  of  wind  tunnel  models)  which  can  usually  be 
avoided  by  the  proper  use  of  turbulence  stimulation [22jt 

Some  instances  have  been  encountered  where  the  prediction  method 
failed  to  converge.  For  example,  in  the  case  of  a body  with  a hemis- 
pherical nose  followed  by  a length  of  parallel  centrebody,  the  method 
only  converged  over  a limited  range  of  transition  position.  In  practice, 
laminar  separation  was  known  to  exist  on  this  body  shape  at  low  Reynolds 
numbers.  The  method  also  failed  to  converge  for  some  bodies  with  a 
length  of  parallel  centrebody  greater  than  about-half  the  total  body 
length  at  Reynolds  numbers  R^  less  than  about  10  . 

Normally  the  computations  were  performed  with  a drag  iteration 
tolerance  of  0.001.  However,  this  had  to  be  increased  to  0.003  for  a 
body  with  a fineness  ratio  VD  = 4 at  Reynolds  number  R^  less  than  10' 
for  reasons  of  stability. 
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The  prediction  method  has  been  modified  to  allow  for  the  effect  of 
a simple,  lightly-loaded  propulsor  and  this  will  be  described  in  a 
future  report.  In  practice,  control  surfaces  are  also  fitted  to  the 
tail  region  which  make  the  flow  highly  three-dimensional.  The  effect 
of  these  cannot  be  incorporated  into  an  axisymmetric  method,  but  the 
additional  drag  resulting  from  these  and  from  surface  roughness  can  be 
estimated  £23}  and  added  to  the  bare-hull  predicted  value. 

2.7  Sources  of  error  in  the  prediction  method. 


This  topic  has  been  discussed  by  Myring  £127,  who  showed  that  any 
error  in  the  level  of  skin  friction  is  likely  to  lead  to  an  error  in  the 
drag  estimation  of  the  same  magnitude.  In  addition  because  of  the 
assumed  relationship  between  skin  friction  and  shape  parameter,  as 
defined  in  referenceflij  a 1%  error  in  shape  parameter  over  the  entire 
body  can  lead  to  a 1.6%  error  in  the  value  of  profile  drag. 

Unfortunately,  accurate  boundary-layer  measurements  on  a body  of 
revolution  are  scarce.  Probably  the  most  detailed  results  available  are 
those  of  Patel  et  al  [19J  obtained  using  a hot-wire  anemometer.  Over  the 
latter  half  of  the  body,  predictions  and  measurements  of  shape  factor  £37 
are  in  good  agreement  except  over  the  last  10%  of  the  body  length.  This 
disparity  will  have  only  a small  effect  on  the  profile  drag.  Hence  any 
error  in  the  prediction  of  profile  drag  is  likely  to  be  due  to  an  error 
in  the  level  of  skin-friction. 

3.  COMPARISON  OF  DRAG  PREDICTIONS WITH  EXPERIMENT 
3.1  Myring 's  comparisons. 

Myring  £l2j  has  made  comparisons  of  drag  predictions  and  measurements 
for  the  airship  'Akron'  £24J  and  two  bodies,  A and  B,  of  Lyon  f25j»  Three 
different^ zed  models  of  'Akron'  gave  different  values  of  drag,  those  of 
the  two  larger  models  being  7%  and  14%  less  than  that  of  the  smallest. 

Both  Young's  and  Myring's  methods  gave  almost  identical  results  which  are 
slightly  larger  than  that  for  the  medium  sized  model.  For  Lyon's  A and  B, 
Young's  method  gave  somewhat  better  agreement  with  experiment,  Myring's 
predictions  being  about  8%  less  than  the  measurements,  tyring  pointed  out 
however  that  Lyon's  measurements  may  be  in  error  by  as  much  as  10%  because 
of  tunnel  blockage  and  other  effects. 

3«2  Comparison  of  predicted  and  measured  values  of  profile  drag. 

The  experimental  data  discussed  in  this  section  were  all  obtained  in 
a situation  where  there  was  no  laminar  separation  and  where  either  the 
location  of  natural  transition  was  known  or  transition  was  fixed 
artificially  by  means  of  turbulence  stimulation.  Predictions  were 
obtained  with  transition  specified  at  the  corresponding  location.  All 
the  predictions  in  this  section  were  made  using  the  Ludwi eg-Ti liman  skin- 
friction  law  £157*  Details  of  the  body  shapes  and  flow  conditions  are 
listed  in  Table  1. 

Predictions  over  a range  of  Reynolds  number  were  obtained  for  three 
streamline  bodies  with  different  fineness  ratios,  referred  to  as  DTMB4154, 
4159  and  4l65,  and  tested  in  a towing  tank  at  David  Taylor  Model  Basin 
(DTMB)  £*267.  The  drag  of  body  4l65  was  also  measured  by  ARL  (now  AMTE 
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(Teddington) ) £"2 7j  mainly  in  low  speed  wind  tunnels,  from  both  drag 
balance  measurements  and  pitot-static  wake  surveys.  In  all  the  measure- 
ments, transition  was  fixed  by  means  of  a trip  wire  or  other  similar 
mechanism.  Comparison  of  predictions  and  measurements  is  illustrated  as 
Figure  l(a).  Balance  and  wake  measurements  were  also  carried  out  on  a 
torpedo-shaped  body  with  a hemispherical  nose  [271  and  these  are  compared 
with  predictions  in  Figure  l(b).  The  agreement  in  both  cases  is  very 
favourable. 

The  effect  of  nose  shape  on  drag  is  shown  in  Figure  2.  Nose  shape 
was  varied  by  changing  the  nose  fineness  ratio  while  keeping  the  nose 
shape  constant  and  the  total  volume  was  kept  constant  by  varying  the 
length  of  parallel  centrebody.  The  profile  drag  of  the  bodies,  referred 
to  as  4620-1  to  4,  was  measured  over  a range  of  Reynolds  number  by 
McCarthy  et  al  [22") • The  location  of  natural  transition  was  determined 
using  hot  film  probes.  Transition  on  the  two  fullest  nose.&bapes  was  a 
result  of  laminar  separation  which  occurred  at  the  same  location  at  all 
the  test  Reynolds  numbers,  but  the  separation  bubbles  so  formed  were 
short  in  length  and  the  parasitic  drag  due  to  separation  was  assumed  to 
be  negligible.  On  the  other  two  bodies,  the  position  of  natural  transition 
moved  forward  with  increasing  Reynolds  number.  The  measurements  of  profile 
drag  are  compared  with  predictions  using  the  present  method,  and  with 
predictions  [22l  using  the  finite  difference  method  of  Cebici  & Smith  [7j> 
The  two  theoretical  methods  give  almost  identical  results.  In  both  cases 
the  trends  are  well  predicted  but  the  actual  values  slightly  under- 
predicted. 

Comparisons  illustrating  the  effect  of  transition  position  on  the 
drag  of  Lyon's  bodies  A and  B£257  (both  streamline  shapes)  is  shown  in 
Figure  3*  Transition  was  controlled  by  varying  the  level  in  the  oncoming 
flow  using  the  presence  or  absence  of  a wire  screen.  The  trends  are 
correctly  predicted  although  the  absolute  values  are  slightly  under- 
predicted. As  mentioned  above  (section  3.1)  the  measurements  may  be 
subject  to  some  experimental  error. 

Correlations  between  predicted  and  measured  values  of  the  profile 
drag  coefficients  C^  (based  on  maximum  cross  sectional  area)  and  Cy 

(based  on  ( volume )^/3)  are  presented  in  Figures  4 and  5«  In  general  the 
agreement  is  as  good  as  the  accuracy  to  which  drag  can  be  measured. 

The  results  in  Figures  4 and  3 are  in  agreement  with  some  results 
presented  by  Chappell  fl67,  who  compared  predictions  using  Myring's  method 
with  the  ESDU  data  sheets  £27.  These  data  sheets  give  C^  as  a function  of 
Reynolds  number  and  fineness  ratio  VD  and  they  were  derived  using  the 
method  of  Young  [ij . Chappell's  conclusion  was  that  in  the  Reynolds  number 
range  2 x 10”<Rjj<2  x 10?  Myring's  method  with  the  Ludwitg-Tillman  skin- 
friction  law  is  likely  to  give  values  of  drag  correct  to  about  - 5$,  a 
reduction  of  about  50$  on  the  likely  errors  using  the  ESDU  data  sheets. 

3.3  Comparison  of  the  Ludwitg-Tillman  and  the  Chappell  skin-friction  laws. 

Figure  6 illustrates  a comparison  between  predictions  of  skin-friction 
drag  coefficient  (based  on  wetted  area)  and  two  correlations  of  flat  plate 
experimental  data.  The  correlations  are  the  ITTC  line  given  by 
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(7) 


Cp  = 0.075  / (logl0BL-2)2 

and  the  Schoenherr  formula  given  by 

CF  = 0.242  / log^HjCp)  2 (8) 

At  low  Reynolds  numbers  the  Schoenherr  formula  deviates  slightly  from 
the  ITTC  line  and  this  is  discussed  by  Granville  f 2817 . The  predictions 
in  Figure  6 are  for  a body  with  a fineness  ratio  (L/D)  of  ten-to-one 
(DTMB  4159)  [ 26 J,  using  Myring's  method  with  both  the  Ludwieg-Tillman 
and  the  Chappell  skin-friction  laws.  Some  computations  with  both  laws 
were  also  made  on  bodies  with  fineness  ratios  of  seven-  and  four-to-one 
(DTMB  4157  and  4154  respectively)  and  the  results  for  the  former  body 
were  close  to  those  of  the  ten-to-one  body.  In  all  cases  transition  was 
prescribed  to  occur  well-forward  (at  x/L  = 0.05).  On  simple  theoretical 
grounds  one  would  expect  the  results  with  increasing  VD  to  approach  from 
above  those  of  a flat  plate  and  these  trends  are  observed  using  the 
Chappell  skin-friction  law.  With  the  Ludweig-Tillman  law  the  correspond- 
ing values  are  10-15#  lower  in  the  Reynolds  number  range  4 x 10°<Rl<109. 
In  some  cases  the  predicted  values  of  profile  drag  coefficient  were  less 
than  the  corresponding  flat  plate  skin-friction  values.  At  Reynolds 
numbers  below  2 x 10?  both  skin-friction  laws  gave  almost  identical 
results. 


4.  COMPARISONS  OF  VELOCITY  PREDICTIONS  AND  MEASUREMENTS 
4.1  Patel  et  al  measurements. 

This  section  contains  comparison  of  the  unmodified  and  modified 
(section  2.4)  prediction  methods  with  the  hot  wire  measurements  of  Patel 
et  al.  C19J  The  test  body  had  a conical  tail  of  about  23°  semi -angle, 
nevertheless  the  flow  remained  attached  right  up  to  the  tip  at  the  test 
Reynolds  number  of  1.26  x 10°.  The  variation  of  velocity  components 
across  the  boundary  layer  were  obtained  at  a number  of  different  axial 
stations  and  these  can  be  used  to  obtain  integral  parameters  which  can 
be  directly  compared  with  the  predictions.  Values  of  skin-friction 
coefficient  were  measured  using  three  different  methods  and  there  was 
reasonable  agreement  between  them. 

Comparison  of  predictions  and  measurements  of  various  integral 
parameters  using  the  Myring  prediction  method  with  and  without  the  blunt 
after-body  modification  mentioned  in  section  2.4  and  equation  6 are  pres- 
ented as  Figures  7 and  8.  Figures  7(a)  and  7(b)  taken  from  reference  3, 
using  the  unmodified  theory,  show  that  the  shape  parameter  H is  under- 
predicted  near  the  tip  of  the  body  although  the  momentum  area  8 is  well 
predicted.  Incorporation  of  the  modification  (Figures  8(a)  and  8(b)) 
leads  to  an  improvement  in  the  prediction  of  shape  parameter  and  skin- 
friction  coefficient  although  the  external  velocity  distribution  and  the 
momentum  area  are  slightly  less  well  predicted. 

Figure  9 shows  a comparison  of  measured  and  predicted  profiles  of 
velocity  component  along  the  hull.  The  improved  shape  parameter 
prediction  clearly  leads  to  an  improved  prediction  of  the  general  shape 
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of  the  velocity  profiles  near  the  tip  of  the  body.  In  particular, 
values  of  velocity  close  to  the  hull  are  better  predicted.  The  agree- 
ment is  good  considering  the  limitations  imposed  by  assuming  a one- 
parameter  family  of  velocity  profiles. 

4.2  Water-tunnel  measurements. 

A number  of  velocity  measurements  have  been  made  on  torpedo-shaped 
bodies  of  revolution  in  the  thirty-inch  water  tunnel  at  AMTE  (Teddington). 
The  bodies  referred  to  in  this  section  all  had  conical  tails  of  semi-angle 
1o  and  flat  noses.  Hemispherical  or  elliptic  noses  were  used  in  the 
boundary-layer  calculations  to  obtain  the  velocity  predictions.  However 
the  effect  of  this  on  boundary-layer  predictions  near  the  tip  of  the  body 
should  be  small. 

Some  velocity  measurements  were  carried  out  on  an  unpowered  torpedo 
model  using  three  rakes  of  total  head  tubes  with  the  static  pressure 
determined  from  static  holes  in  the  hull.  The  measurements  were  made  in 
planes  normal  to  the  hull  corresponding  to  those  of  the  predictions  at 
two  different  Reynolds  numbers.  Comparison  of  the  measured  total  velocity 
profiles  with  the  predicted  velocity  components  along  the  hull  is  shown 
in  Figure  10.  There  is  some  scatter  on  the  measurements  and  the  agree- 
ment between  measurements  and  predictions  deteriorates  with  downstream 
position.  One  possible  reason  for  this  disparity  is  the  error  introduced 
in  the  measurements  by  assuming  that  the  static  pressure  does  not  change 
across  the  boundary  layer.  As  one  moves  downstream  there  may  be  a 
variation  in  static  pressure  across  the  boundary  layer  which  would 
introduce  an  error  in  the  measurements . Another  possible  reason  is  that 
measurements  of  total  velocity  are  being  compared  with  predictions  of 
velocity  component  along  the  hull.  The  greatest  discrepancy  is  likely 
to  be  at  the  edge  of  the  boundary  layer  where  the  flow  is  axial,  in 
which  case  the  predictions  ought  to  be  about  4#  less  than  the  measurements. 
At  the  most  downstream  position  the  disparity  is  in  excess  of  10#  across 
the  whole  boundary  layer,  which  means  that  the  latter  correction  is  not 
sufficient  to  account  for  it. 

A single  boundary  layer  traverse  was  carried  out  on  a torpedo-shaped 
body  similar  to  the  previous  one,  but  fitted  with  a fin-duct  assembly. 
Velocity  was  measured  using  a pitotstatic  tube  a short  distance  upstream 
of  the  duct  leading  edge  and  in  a plane  normal  to  the  axis.  Predictions  in 
this  plane  were  obtained  by  assuming  that  the  velocity  component  parallel 
to  the  hull  does  not  change  along  a line  drawn  parallel  to  the  hull  between 
the  plane  of  the  predictions  and  that  of  the  measurements.  From  Figure  11 
in  which  predictions  and  measurements  are  compared,  it  is  evident  that 
there  is  little  scatter  on  the  measurements.  The  predictions  labelled  as 
'uncorrected ' in  Figure  11  agree  with  the  measurements  near  to  the  hull, 
but  they  deviate  slightly  with  increasing  distance  from  hull.  This 
disparity  can  be  largely  explained  by  considering  that  the  measurements 
are  of  total  velocity  and  predictions  are  of  velocity  component  along  the 
hull.  Near  the  hull  the  streamlines  will  be  parallel  to  the  hull  but  at 
the  edge  of  the  boundary  layer  the  total  velocity  will  be  almost  axial  in 
direction.  Between  these  extremes  the  direction  of  flow  will  vary  between 
the  two  values.  To  obtain  the  total  head  value,  the  component  along  the 
hull  should  be  divided  by  the  cosine  of  the  angle  between  the  flow  and  the 
hull  and  vice  versa.  Figure  11  shows  predicted  values  of  total  velocity 
calculated  this  way,  assuming  a uniform  variation  of  flow  angle  between 
the  hull  and  the  predicted  edge  of  the  boundary  layer.  The  agreement  is 
now  improved.  At  the  edge  of  the  boundary  layer,  a flow  angle  of  l6°, 
the  cgrrection  to  the  predicted  value  is  about  4#.  At  a corresponding  angle 
of  12°  the  correction  would  be  about  2#  and  at  20°  about  6#. 
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4.3  Wind-tunnel  measurements. 

A number  of  wind  tunnel  measurements  have  been  made  on  bodies  of 
revolution  by  AMTE  (Teddington)  in  the  NMI  low-speed  wind  tunnels.  Two 
different  bodies  were  tested,  one  a streamline  body  with  a maximum  tail- 
cone semi-angle  of  about  13°  and  the  other  having  a length  of  parallel 
centrebody  with  a much  fuller  conical  after-end  with  a semi-angle  of 
about  22°.  The  hull  static  pressure  distribution  was  measured  and  wake 
surveys  were  conducted  at  one  axial  station  near  the  tip  on  the  former 
body  and  two  on  the  latter,  using  a pitotstatic  probe  connected  to 
traversing  gear  programmed  to  move  in  radial  and  circumferential 
directions.  These  surveys  showed  the  flow  to  be  axisymmetric  apart 
from  the  region  downstream  of  a thin  supporting  strut. 

Figure  12  (taken  from  reference  3)  shows  comparisons  between  the 
measured  and  predicted  external  velocity  distributions,  which  are 
related  to  the  corresponding  pressure  distributions.  Figure  13  shows 
comparison  between  velocity  measurements  and  predictions  at  one  axial 
position  on  the  streamline  body,  corrected  in  the  manner  described  in 
section  4.2  so  as  to  be  of  the  total  velocity  along  a radius.  Although 
the  shape  parameter  is  slightly  under  predicted  as  mentioned  in  reference 
3,  the  agreement  between  values  of  velocity  is  good. 

A further  series  of  measurements  was  carried  out  on  the  full-after- 
• end  body.  Four  rakes  which  could  be  either  total  or  static  head  ..tubes 

were  fixed  to  the  cailcone  at  different  stftions.  The  two  middle 
stations  corresponded  to  those  of  the  programmed  pitotstatic  probes 
. mentioned  above.  Figure  14  illustrates  comparison  between  predictions 

and  measurements  at  all  four  rake  positions.  Predictions  were  obtained 
with  and  without  the  modification  described  in  section  2.4.  The  measure- 
ments using  the  total  and  static  head  tubes  show  considerably  more 
scatter  than  those  using  the  programmed  pitotstatic  tubes.  At  the  second 
station  the  latter  measurements  agree  slightly  better  with  the 
predictions  obtained  with  modification.  At  the  next  station  downstream 
neither  prediction  agrees  well,  but  the  unmodified  method  is  in  slightly 
better  agreement.  The  measurements  clearly  indicate  that  the  flow  has 
separated  at  the  most  downstream  position. 

Some  recent  measurements  were  conducted  on  the  body  with  full  after- 
end with  a lightly  longer  length  of  parallel  middle  body  in  the  AMTE 
(Haslar)  Ho  2 Ship  Tank.  A rake  with  twelve  pitotstatic  tubes  was  used 
to  measure  total  velocity  at  x/L  = 0.985  corresponding  to  the  measuring 
position  furthest  downstream  in  the  wind  tunnel  tests,  but  at  about  twice 
the  Reynolds  number.  At  this  higher  Reynolds  number  no  separation  is 
present,  but  the  unmodified  theory  (Figure  15)  gives  slightly  better 
agreement  with  the  measured  velocity .profiles.  Figures  14  and  15  illus- 
trate that  near  the  tip  of  a full  afterbody,  where  the  flow  has  either 
already  separated  or  is  close  to  separation,  the  power  law  assumption 
(equation  4)  breaks  down.  Also,  broadly  speaking,  the  modified  theory 
gives  better  velocity  predictions  near  the  hull  and  the  unmodified  theory 
gives  better  predictions  near  the  edge  of  the  boundary  layer. 

In  practice,  the  tailcone  of  a body  of  revolution  is  often  truncated. 
This  will  cause  separation  of  the  flow  if  it  has  now  already  occurred  which 
may  well  have  an  effect  on  the  flow  for  some  distance  upstream.  The 
prediction  method  assumes  that  the  tailcone  ends  in  a point.  The  difference 
may  account  for  some  of  the  above  disparities. 
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5.  CONCLUSIONS  AMD  SUGGESTIONS  FOR  FURTHER  RESEARCH 


In  general,  the  drag  prediction  method  of  Myring  with  the  akin-friction 
law  of  Chappell  givea  predictions  which  agree  well  with  measurements  both  in 
trend  and  magnitude.  The  predicted  trends  correspond  with  those  of  other 
theoretical  and  semi -empirical  methods.  The  method  has  now  been  used  by  ESDU 
to  produce  data  item  7 8019  £*29J  on  the  profile  drag  of  axisymmetric  bodies 
at  zero  incidence  for  subcritical  Mach  numbers. 

In  general,  velocity  measurements  and  predictions  using  the  unmodified 
Myring  method  agree  moderately  well,  except  near  the  tip  of  a full  after-body. 
In  this  region,  especially  at  low  Reynolds  numbers,  the  flow  is  often 
separated  or  close  to  separation  and  the  assumption  of  a power  law  velocity 
profile,  inherent  in  the  method,  breaks  down.  The  modified  theory  of  section 
2.4  predicts  values  of  velocity  near  the  hull  rather  better  than  the 
unmodified  theory,  but  the  position  is  reversed  at  the  edge  of  the  boundary 
layer.  In  general,  the  unmodified  theory  gives  better  velocity  predictions 
than  the  modified  theory.  It  is  hoped  to  be  able  to  use  the  method  to  scale 
model  test  data  on  a fully  appended  body  to  full  scale,  where  the  difference 
in  Reynolds  numbers  may  be  a factor  of  as  much  as  a hundred. 

At  present  the  prediction  method  breaks  down  on  bodies  with  flat  noses, 
because  the  laminar  flow  method  becomes  unstable  in  large  favourable  or 
adverse  pressure  gradients.  However,  there  is  experimental  evidence  which 
shows  that  these  are  bodies  with  flat  noses  where  laminar  separation  does  not 
occur  over  some  range  of  Reynolds  number.  The  method  could  be  modified  to 
deal  with  a body  with  a flat  nose  by  starting  the  laminar  boundary-layer 
calculation  at  the  end  of  the  nose  flat  using  boundary-layer  parameters 
derived  from  a separate  calculation  using  a differential  laminar  boundary- 
layer  method. 

In  any  real  situation  a body  of  revolution  will  normally  be  fitted  with 
control  and  stabilizing  fins  and  a propulsor.  The  Myring  prediction  method 
has  been  modified  to  include  the  effect  of  a single  propeller  on  the  boundary- 
layer  development  and  this  modification  will  be  described  in  a future  report. 
There  are  empirical  methods  available  for  estimating  the  drag  increment  due  to 
fins  and  the  thrust  deduction  of  the  propulsor,  but  the  accuracy  of  these 
estimates  and  how  they  scale  is  subject  to  considerable  uncertainty.  Neverthe- 
less the  present  prediction  method  should  prove  useful  in  feasibility  and 
design  studies. 
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FIG.  4 CORRELATION  BETWEEN  MEASURED  AND  PREDICTED  VALUES  OF 
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FIG.  5 CORRELATION  BETWEEN  MEASURED  AND  PREDICTED  VALUES  OF 
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FIG. 6 COMPARISON  OF  PREDICTIONS  USING  DIFFERENT 
SKIN  FRICTION  LAWS  WITH  KNOWN  CORRELATIONS 
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FIG.  7(b)  COMPARISON  OF  PREDICTIONS  OF  INTEGRAL  PARAMETERS 
WITH  THE  MEASUREMENTS  OF  PATEL.  ET  AL  119)  (UNMODIFIED  THEORY) 
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PREDICTIONS 


VELOCITY  PREDICTIONS  ANO  MEASUREMENTS  ON  PATEL  ET  AL  BODY  [19] 


VELOCITY  PREDICTIONS  AND  MEASUREMENTS  ON  UNPOWERED  TORPEDO  MODEL 
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VELOCITY  PREDICTIONS  AND  MEASUREMENTS  ON  A 
STREAMLINE  800Y  OF  REVOLUTION 
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FIG.  15  VELOCITY  PREDICTIONS  AND  MEASUREMENTS  ON  A 
BODY  OF  REVOLUTION  WITH  FULL  AFTER- END 
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